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We report the experimental observation of novel vortex patterns in vertically vibrated granular 
rods. Above a critical packing fraction, moving ordered domains of nearly vertical rods sponta- 
neously form and coexist with horizontal rods. The domains of vertical rods coarsen in time to form 
large vortices. We investigate the conditions under which the vortices occur by varying the number 
of rods, vibration amplitude and frequency. The size of the vortices increases with the number of 
rods. We characterize the growth of the ordered domains by measuring the area fraction of the or- 
dered regions as a function of time. A void filling model is presented to describe the nucleation and 
growth of the vertical domains. We track the ends of the vertical rods and obtain the velocity fields 
of the vortices. The rotation speed of the rods is observed to depend on the vibration velocity of the 
container and on the packing. To investigate the impact of the direction of driving on the observed 
phenomena, we performed experiments with the container vibrated horizontally. Although vertical 
domains form, vortices are not observed. We therefore argue that the motion is generated due to 
the interaction of the inclination of the rods with the bottom of a vertically vibrated container. We 
also perform simple experiments with a single row of rods in an annulus. These experiments directly 
demonstrate that the rod motion is generated when the rods are inclined from the vertical, and is 
always in the direction of the inclination. 



I. INTRODUCTION 

Granular materials are well known examples of dissi- 
pative nonequilibrium systems that show a rich variety of 
collective phenomenon such as convection, wave patterns, 
and segregation [1-4]. Most studies utilize spherical par- 
ticles to investigate these bulk properties. However, in 
most natural or industrial settings one can find an abun- 
dance of anisotropic granular materials viz. rice, medicine 
capsules, and even logs. Therefore it is surprising that 
very few studies of prolate granular materials have been 
carried out. Mounfield and Edwards [5] applied the con- 
cepts of configurational statistical mechanics to study the 
nature of the isotropic to nematic phase transition in 
a granular system of elongated particles. In recent ex- 
periments utilizing a tall narrow cylinder, Villarruel et 
al. [6] studied the effects of anisotropy on granular pack- 
ing. They observed the appearance of smectic states with 
the direction given by the container walls. 

In thermal systems, particle anisotropy is known to 
produce ordered states. Examples are rod and plate 
shaped colloids and liquid crystals, which show orient a- 
tional order and form nematic and smectic phases [7-11]. 
The ordering mechanism was found to be entropically 
driven, i.e. as thermodynamic equilibrium is reached, the 
process of entropy maximization leads to long range or- 
der. It is not obvious that such mechanisms carry over 
to granular systems because the thermal energy scale is 
very much smaller than the potential energy required for 
particle rearrangement, and energy has to be supplied 
actively to produce sustained motion. Therefore, an in- 
teresting question arises - does shape anisotropy lead to 
self-organization and pattern formation in granular ma- 
terials? 

In this paper, we report the observation of novel vor- 
tex patterns exhibited by granular rods that are vibrated 



vertically inside a container. We obtain the phase dia- 
gram for the observed patterns as a function of the ac- 
celeration of the driving and the packing fraction of the 
rods. We find that for sufficiently large packing fractions, 
the rods tend to align vertically and undergo vortex mo- 
tion. Using high frame rate imaging and particle track- 
ing, we have measured the velocity fields of the vortices 
as a function of packing fraction and driving frequency. 
Based on our experimental observations, we argue that 
the inclination of the rods causes motion due to collisions 
with the bottom boundary. To bolster our claim, we con- 
duct experiments with single row of rods in an annulus 
to demonstrate that motion is always in the direction of 
inclination. Thus we show that shape anisotropy leads 
to translational motion in such systems. 



II. EXPERIMENTAL APPARATUS 

The experiments were performed in a circular anodized 
aluminum container with a diameter D = 6.0 cm and 
depth H = 1.5 cm. (Limited experiments were also car- 
ried out in a container with D = 8.5 cm.) The con- 
tainer is leveled to within 0.002 cm and is attached to an 
electro-mechanical shaker through a rigid linear bearing 
that allows motion only in the vertical direction. The 
cell is driven by a sinusoidal signal at a frequency /, and 
is monitored by an accelerometer via a lock-in amplifier. 
The experiments were conducted using copper cylinders 
with uniform length / = 6.2 mm and diameter d = 0.5 
mm. The cylindrical surface of the rods is coated with a 
grey tin-oxide layer to diminish light reflection and has 
a coefficient of friction /i = 0.32. The patterns that form 
are imaged from above using a high-frame rate digital 
camera (Kodak SR-1000). Since the flat ends of the rods 
reflect light better than the sides, they appear as bright 
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FIG. 1: The phase diagram of the observed patterns in vi- 
brated granular rods. The driving frequency / = 50 Hz. Vor- 
tices are observed for sufficiently large numbers of rods and 
driving acceleration parameter T. 



spots when imaged from above. If the rods are inclined 
greater than 35 degrees, they reflect far less light com- 
pared to nearly vertical and horizontal rods. 

One control parameter for our experiments is T = 
A(27rf) 2 /g, where A is the driving amplitude and g is the 
acceleration due to gravity. A second control parameter 
is the non-dimensional number fraction = N/N maxi 
where N is the total ^ number of rods in the container 
and N max = -J= is the maximum number of rods 

required to obtain a vertically aligned monolayer. There- 
fore, = 1 corresponds to having one layer of triangu- 
larly packed vertical rods. 




FIG. 2: Examples of patterns observed when granular rods 
are vertically vibrated (/ = 50 Hz), (a) Nematic- like state 
(n^ = 0.152, r = 2.42). (b) Moving domains of near vertical 
rods (n0 = 0.344, Y = 3.38). (c) Multiple vortices {n<f> = 
0.551, T = 3.29). (d) Large vortex (n = 0.535, F = 3.00). 
The cylindrical surface of the rods is coated with a grey tin- 
oxide layer to diminish light reflection and has a coefficient 
of friction \i — 0.32. The patterns that form are imaged 
from above using a high-frame rate digital camera (Kodak SR- 
1000). Since the tips of the rods reflect light better than the 
sides, they appear as bright spots when imaged from above. 
If the rods are inclined greater than 35 degrees, they reflect 
far less light compared to nearly vertical and horizontal rods. 



III. OBSERVED PATTERNS 

The system is initialized by pouring the rods into the 
container and then increasing T, which establishes a ran- 
dom state. Figure 1 shows the phase diagram of the 
observed patterns as a function of (f = 50 Hz). The 
rods are observed to remain static in the initial configu- 
ration until T > 1.47. As T is increased above this value, 
the rods heap towards one side of the container similar to 
previous observations with spherical particles [12] . As the 
acceleration is increased further above T = 2.6, the rods 
are observed to spread-out evenly inside the container, 
and form horizontal layers with random orientations. An 
example of such a nematic phase observed at low n$ is 
shown in Fig. 2(a). Above a critical n^, a second transi- 
tion is observed. Chaotically moving domains of almost 
vertical rods are observed to spontaneously form and co- 
exist with horizontal rods [13]. An example is shown in 
Fig. 2(b). 

As is increased further, the domains of near- vertical 
rods that are formed coalesce and undergo vortex motion. 
Fig. 2(c) shows an example of two counter-rotating vor- 
tices. For large number fractions (n<j, > 0.49), a single 



large stable vortex made entirely of near- vertical rods in 
the center and inclined rods at its' boundary, is observed 
[Fig. 2(d)]. 

If the acceleration is increased further, a gas like state 
is reached where the rods vibrate vigorously inside the 
container and the nematic and vortex states are de- 
stroyed. For the highest n^, a pure gaseous state is not 
reached due to the limits of the apparatus. Instead, do- 
mains of ordered vertical rods are observed to coexist 
with the gas of rods. 



A. Domain Growth 

We now discuss how vortices nucleate and grow as a 
function of time. An example of the growth process is 
shown in Fig. 3, and a movie can be viewed at Ref. [14]. 
Starting from a random state, pockets of near-vertical 
domains are observed to nucleate uniformly inside the 
container. This is in contrast with previous observations 
where the ordered domains grow inward from the bound- 
ary [6]. The pockets of vertical domains grow in time by 
merging to form larger domains of almost vertical rods. 
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FIG. 3: The formation of the vortex from an initial random 
state as function of time t of vibration with r = 2.6 and 
ricf) — 0.39. (a) t = 3 min, (b) t = 6 min, (c) t = 9 min, (d) 
t = 12 min. Small domains of vertical rods are observed to 
form and coarsen over time and evolve into a vortex. (/ = 60 
Hz). 

The domains are then observed to collectively move as 
the the system becomes more ordered and then finally 
show vorticity. 

We measure the growth of the domains by taking time- 
lapse images of the process described above. The growth 
of near- vertical domains is then measured from the ratio 
of the domain area to that of the total container area as 
a function of time. The fraction of rods that are near- 
vertical, n v is plotted in Fig. 4 for a range of at a fixed 
driving acceleration and frequency, T = 3.4 and / = 50 
Hz. 

We find that there exist two growth regimes [a model 
is discussed in Sec. V]. For short times, the growth rate 
depends on in the following way. If is below a 
critical value, domains appear but never organize into 
large vortices. However, as is increased the short time 
growth is slowed by the high packing. That is, more time 
is required for small domains to form due to the lack of 
voids present. A void filling model for domain growth 
will be discussed in Sec. V. The saturation value of n v 
increases with n^. The difference in the values arises due 
to the definitions of n v and n^. 

The variation of the relative number of near-vertical 
rods as a function of is plotted in Fig. 5. We show 
separate data for the relative number of inclined rods, 
vertical rods and horizontal rods. Rods that are tilted 
between 30 and 80 degrees are considered inclined, and 
the rest are considered as either horizontal or vertical 
(the inclination of the rod was obtained by measuring 
the distance x between adjacent rods and calculating 
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FIG. 4: The fraction of vertical rods Tin as a function of time 
for various at T = 3.4 and / = 50 Hz. The vertical do- 
mains show two distinct growth regimes. At short times, the 
growth is increasing exponentially and then gives way to an 
exponential decrease in growth to a saturation value deter- 
mined by n^. For — 0.38, the domain growth saturates 
very slowly to a value less than because the near verti- 
cal domains cannot be supported by the horizontal rods. For 
> 0.46 the two distinct growth regimes are apparent and 
for all 7i(f) vortex motion occurs after t = 300 s. We note 
that for the highest value of n<j> the initial growth displays a 
slow down, which is directly related to the void filling mech- 
anism discussed in the text. The dashed line is a fit to the 
integrated growth equation discussed in Sec. V. We find that 
the simplified model gives an accurate interpretation of the 
results. 

sm~ 1 [d/x]). For < 0.22, the container has only hor- 
izontal rods. As increases, more rods tend to the 
vertical direction. For > 0.71 most rods are either 
vertical or horizontal. 



B. Spatial structure of the vortex and velocity 
fields 

We next discuss the spatial structure of the vortex. In 
Fig. 6, a close-up image of a vortex pattern is displayed 
which shows the progressive inclination of the rods from 
the center of the vortex. We note that the direction of 
motion corresponds to the direction of inclination (for 
Fig. 6 the motion is counterclockwise). In this case the 
inclination of the rods varies from approximately 85 to 
40 degrees near the edges. The change in inclination of 
the rods within the vortex is less at higher due to the 
increased packing. 

By acquiring images at 250 (frames s _1 ), we were able 
to track the ends of the vertical rods and obtain the ve- 
locity fields of the vortices. In addition to the rotational 
motion, rods also vibrated due to the collisonal interac- 
tions with their neighbors. We first track [15] the ends 
of the near-vertical rods and then perform spatial and 
temporal averaging to obtain the vortex fields. An ex- 
ample of the velocity field of a vortex is presented in 
Fig. 7 (here the temporal averages include 100 consecu- 
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FIG. 5: The relative number of inclined, vertical and hori- 
zontal rods as a function of n$ after steady state is reached 
r 3.00 and / = 50 Hz. 



tive frames, and spatial averages were taken over an area 
of 2d x 2d). 

The averaged azimuthal velocity v(r\ as a function of 
the distance r from the center of the vortex, is shown in 
Fig. 8(a) for a range of n^. The data was obtained at a 
constant acceleration T = 3.00 for a frequency / = 50 Hz. 
We note that for small distances relative to the center of 
the vortices, the averaged velocity v(r) increases linearly 
with the distance r indicating solid body rotation. At 
intermediate ranges, this linear relationship is not ob- 
served, and therefore shearing occurs inside the vortices. 
As the boundary of the vortex is approached, velocity 
decreases due to friction with the horizontal rods at the 
edge. It is interesting to note that the horizontal rods 
at the edge of the vortex are often aligned tangentially 
with the boundary of the vortex [see Fig. 2(c)]. We also 
observe that the slope of the averaged velocity v(r) at 
small r systematically decreases as the packing fraction 
is increased. Thus the angular velocity of the inner 
core of the vortex decreases with its size. 



C. Frequency dependence 

We also explored the rotation rate dependence of the 
vortices as a function of /. Figure 8(b) shows v(r) for 
= 0.53 and T = 3.00 for vortices of equivalent size. 
As stated, F is held constant while frequency is increased, 
which decreases vibration velocity. Therefore, we demon- 
strate that the vortex speed systematically increases with 
vibration velocity. 

IV. HORIZONTALLY VIBRATED CONTAINER 

To test if the transition to a vertical state is dependent 
on the way the container is vibrated, experiments were 
performed with horizontal driving. There do in fact exist 
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FIG. 6: Close-up image of a vortex that has self organized 
at rid) — 0-419 (only the upper half of the vortex is shown). 
The rods are inclined in the direction of rotation, which is 
counterclockwise. The inclination angle of the rods varies 
from about 90 degrees at the center of the vortex to near 40 
degrees close to the edge. 



similarities: at low heaping and nematic like domains 
of horizontal rods are observed; at high the rods tend 
to align vertically and form ordered domains. An image 
is shown in Fig. 9. In contrast to the vertical shaking the 
vertically aligned domains do not migrate and coarsen 
into vortices. We also find that in this case the convec- 
tion is stronger and destroys the motion of the vertically 
ordered domains. 



V. DISCUSSION 

We first discuss possible mechanisms responsible for 
the rod to align vertically at high packing fractions. The 
tendency of the rods to align vertically at high packing 
fractions may be understood in terms of a void filling 
mechanism. If voids (i.e. space between rods), are small 
then they can only be filled with a near- vertical rod. On 
the other hand, larger voids will accommodate a horizon- 
tal rod. Assuming that the distribution of voids decreases 
with void size, the most probable configuration will be re- 
gions of vertical rods. Furthermore, the decrease in the 
number of large voids is enhanced by the coalescence of 
regions of vertically aligned rods. This also drives the 
system to pack more closely and produces a decreased 
overall center of mass. Thus ordering at high packing 
fractions can be explained by a process that is analogous 
to configurational entropy driven ordering in thermal sys- 
tems [5]. 

To describe the nucleation and growth of the vertical 
domains at high shown in Fig. 3 and 4 we present 
a simple analysis. The rods are assumed to be either 
vertical or horizontal. We then assume that the growth 
rate is initially linear and asymptotically must decrease 
to zero as the number of horizontal rods is diminished 
and a steady state is reached. Then the evolution of n v 
is described by, 



FIG. 7: The velocity field of a granular vortex (n^ = 
0.497, r = 3.00). The data was obtained by tracking indi- 
vidual rod ends and then by spatial and time averaging the 
obtained displacements (see text). 



= an v (p - n v ). (1) 

where, a and (3 are constants that depend on n^. Eq. 1 
can be integrated and is fit to the measured n v (t) in 
Fig. 4. This simplified interpretation seems to cap- 
ture the nucleation, growth, and saturation of the near- 
vertical domains in our experiments. 

In previous work, Villarruel et al [6], the experiments 
were carried out in a system where the diameter of the 
container and the length of the rods were comparable. 
The container was tapped and thus the rods experienced 
considerable shearing with the side boundary. They ob- 
served that vertical domains nucleated at the boundaries 
and subsequently propagate to the center of the cell. 
These observations are consistent with our findings. We 
have shown that vertical domains can nucleate away from 
the sidewalls and form independent of the direction of the 
driving. 

Next, we elucidate the physical mechanism that is re- 
sponsible for the vortex motion. We performed addi- 
tional experiments with a row of cylindrical rods with 
I = 5.1 cm and d = 0.6 cm in a 1-D annulus with a mean 
radius of 5.5 cm [see Fig. 10]. When subjected to ver- 
tical motion, the rods were always observed to move in 
the direction of the inclination. No translation motion is 
observed when the rods were vertical. A movie showing 
this property can be found at Ref. [14]. We also used ball 
point pens and pencils to check that the detailed shape 
of the rod and the tip is not important to this trans- 
lation mechanism. We find that the translation speed 
depends on the driving frequency and the inclination, 
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FIG. 8: (a) The azimuthal averaged velocity v(r) as a function 
of the distance r from the center of the vortex. The velocity 
is observed to increase linearly from the center of the vortex 
and then reach a maxima before decreasing to the edge of 
the vortex. The rate of increase of velocity at the center is 
observed to systematically decrease with n$ and hence vortex 
size, (b) v(r) as function of frequency / for a fixed T at 
= 0.53. The velocity of rotation is observed to depend on 
the velocity of the vibration. F — 3.00 for all data shown. 



qualitatively similar to that observed for the vortices. 

The physical mechanism for the motion of the inclined 
rods, based on our observation, appears to be as follows. 
When the inclined rods are vibrated vertically, they hit 
the bottom plate at a point away from their center of 
mass. Because the rotation of the rod about its center of 
mass is constrained due to the neighboring rods, it gets 
launched in the direction of the inclination. Thus, the 
greater the inclination of the rods, the further they get 
launched and land (viz. projectile motion), giving rise to 
translational motion. 

From these direct observations we conclude that the in- 
clined rods form the engine that drives the vortex. The 
vertical rods in the center of the vortex are simply pulled 
around by shear induced by the inclined rods. As the 
number of rods is increased and the size of the vortex 
increases, the inclination of the rods also decreases be- 
cause of greater packing. Therefore, the speed of the very 
large vortices is expected to decrease with the size of the 
vortex, consistent with our observation [see Fig. 8(a)]. 
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FIG. 9: Image of domains of vertical rods formed when the 
container is vibrated horizontally. We note that vortices are 
not observed. The direction of vibration is from left to right, 
r = 1.04, / = 30 Hz, and n = 0.612. 




FIG. 10: Image of rods contained in a vertically vibrated 
1-D annulus. The rods are always observed to move in the 
direction of inclination. See text and Ref. [14] for more details. 



We also note that the frequency dependence of v(r) is 
compatible with the mechanism described above. By de- 
creasing the frequency of the driving signal at fixed T 
the velocity of the driving plate is increased. Therefore, 
one expects v(r) to increase with decrease in frequency, 
consistent with our observations [see Fig. 8(b)]. 



VI. CONCLUSION 

In conclusion, granular rods are observed to form ver- 
tically aligned domains at high packing fractions, inde- 
pendent of the direction of vibration. Novel vortex pat- 
terns are observed when the rods are vibrated vertically. 
We have measured the growth of near- vertical domains 
and have developed a simple void filling model that well 
describes our results. We have also shown a new trans- 
lation mechanism which occurs due to anisotropy. Based 
on these observations, Aranson and Tsimring [16] have 
developed a phenomenological model that describes the 
formation and coarsening of the vortices. 
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